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Abstract. The Galactic star forming region in the southern 
sky, associated with IRAS 09002-4732 has been mapped si- 
multaneously in two far infrared bands (A e / / = 148 & 209 /itn), 
with ~ 1' angular resolution. Fifteen sources, including IRAS 
08583-4719, 08589-4714, 09002-4732 and 09014-4736 have 
been detected, some of which are well resolved. Taking advan- 
tage of similar beams in the two bands, a reliable dust tempera- 
ture [T( 148/209)] map has been obtained, which detects colder 
dust (< 30 K) in this region. The HIRES processed IRAS maps 
at 12, 25, 60 and 100 /im, have also been used for comparison. 
The optical depth maps, T200 and T100, generated from these 
far-IR data quantify the spatial distribution of the dust. The dif- 
fuse emission from this entire region has been found to be 35 
% of the total FIR luminosity. The slope of the IMF in the mass 
range 4-16 M© has been estimated to be — 1.25lg'g5 f° r m i s 
star forming complex. 

Radiative transfer models in spherical geometry have been 
explored to fit available observations of the 4 IRAS sources and 
extract various physical parameters for corresponding dust-gas 
clouds. Whereas a constant (r°) radial density distribution is 
favoured in IRAS 08583-4719, 08589-4714 and 09002-4732, 
the r- 1 law is inferred for IRAS 09014-4736. The dust com- 
position is found to be similar (Silicate dominated) in all the 
four IRAS sources modelled. The luminosity per unit mass is 
found to be in the narrow range of 44 - 81 L /M Q for these 
star forming regions. 

Key words: ISM : H II region - ISM : dust - ISM : individual 
: IRAS 09002-4732- Infrared : ISM : continuum 



1. Introduction 

The study of far infrared (FIR) emission from the relatively 
denser regions of the Galactic interstellar medium (ISM) is 
useful for understanding the star formation process. High spa- 
tial resolution mapping of Galactic molecular clouds in the FIR 
continuum (in > 2 bands) provides direct information about the 
physical conditions (e.g. total luminosity, dust distribution, dust 



temperature etc.) of the region. For a constant gas-to-dust ratio, 
the dust optical depth traces matter along the line of sight. A 
comparison of the dust temperature distribution with that of the 
gas temperature (obtained from molecular line measurements) 
can give important information on the dust-gas coupling in star 
forming regions. For medium to high mass embedded YSOs, 
dust is the dominant source of heating of the gas (outside the 
H II region). A FIR study of an individual star formation com- 
plex (at a known distance) can measure the luminosity function 
above a threshold which can usually be translated to the present 
day mass function and even initial mass function. 

Systematic differences can be expected between star form- 
ing regions at different galactocentric distances, since the star 
formation process depends on the physical conditions of the 
local ISM (e.g. density, composition, radiation field etc.). This 
important issue has been addressed by many authors using var- 
ious techniques. Systematic studies in radio continuum and 
molecular line emission based on colour selected samples from 
the IRAS Point Source Catalog (hereafter IRAS PSC) have led 
to very useful databases which can be the starting point for fur- 
ther studies (Wood & Churchwell [19894 Wouterloot & Brand 



1989) 



Although the IRAS mission has proved to be a corner stone 
in improving our understanding of star formation in the Galaxy, 
its sub-optimal angular resolution in the 60 and 100 /im bands 
has been a major limitation. In addition, recent ISO measure- 
ments of nearby spiral galaxies have shown the importance of 
emission at wavelengths > 100 /im and colder dust (T < 25 K) 
seems to be a dom inant constituent of spiral galaxies in general 
(Alton et al. 1998 ). This suggests a similar role for colder dust 
in our Galaxy also. 

In view of the above, the TIFR 1 -metre balloon borne tele- 
scope is flown regularly to map Galactic star forming regions 
simultaneously in two FIR bands with near diffraction limited 
angular resolution. These FIR bands have been chosen to be 
trans-IRAS (A e // = 148 & 209 /mi), so that the distribution 
of interstellar dust cooler than « 30 K which might have been 
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missed by IRAS, can be traced (Ghosh 2000, Mookerjea et al. 
1999[ [20004 |2000b[ ). Accessibility of more than 90 % of the 
Galactic plane (I = 315° - 0° - 290°) from the Indian balloon 
launching station has been exploited by targeting the southern 
sky, which is observationally less explored in the FIR. 
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The catalogue of CO emission associated with colour se- 
lected IRAS PSC sources by Wouterloot & Brand (|1989|; here- 



after WB89) covers the longitude range of the Galactic plane 
of interest to us. Based on information from this catalogue, the 
star forming complex around IRAS 09002-4732 was selected 
for the present study because : (i) all the four bright IRAS PSC 
sources in this region have similar kinematic distances (WB89; 



Sugitani & Ogura 1994); three of them (IRAS 08589-4714, 



09002-4732 and 09014-4736) are known to be physically as- 
sociated, hence this complex presents an opportunity to derive 
its PDMF / IMF if a reasonable number of sources are detected; 
(ii) it lies marginally outside the solar circle, hence a compari- 
son of the physical conditions in this region with that in the 4-5 
kpc molecular region and regions inside the solar circle would 
be interesting; (iii) the complex is very close to the Galactic 
mid plane with z ss 20 pc (WB89) as compared to the scale 
heig ht of 90 pc for the embedded OB stars (Wood & Church- 
well |!989b| ). This region around IRAS 09002-4732 is associ- 
ated with the very extended Vela Molecular Ridge detected in 
the wide angle survey of 12 CO emission by May et al. (1988) 



and Murphy and May (1991), and more recently mapped by 
Yamaguchi et al. ( |1999| ) in both 12 CO and 13 CO lines. 

Most of the available observational data for this region per- 
tain to IRAS 09002-4732 itself and clearly indicate the pres- 
ence of an embedded young stellar object in a dense molec- 
ular cloud. These include measurements of radio continuum 



(Casw ell & Haynes |1987t Manchester & Goss |1969| ; Whiteoak 
1992j; Walsh et al. |1998|), infrared spectra from the IRAS 



Low Resolut ion Sp ectrometer Catalog (hereafter LRS; Olnon 
& Raimond, 1986 ), ISO Long Wavelength Spectrometer data 
(ISO-LWS; Clegg et al. 1996| ) and photometric and i maging 



measurement s in t he near and mid infrared (Lenzen |1991 
Walsh et al. 19991. In addition, millimeter wave continuum 



emission (Reipurth et al. 



1996|) and molecular line detection 
A recent 



1995) 



of CS have been reported (Zinchenko et al 
multi transitional study of IRAS 09002-4732 has led to the de- 
tection of bipolar outflow emission in the CS line with distinct 
red and blue shifted peaks displaced symmetrically relative to 
the central IR object (Lapinov et al. 1998 ). However, no bipolar 
emission in the CO line has been seen by them. 

The observations of a large region (area > 500 sq. arc min) 
centered at IRAS 09002-4732 are presented here. This Galactic 
H II region / molecular cloud complex has been studied with 
the aim of understanding the energetics, physical sizes, spatial 
distribution of interstellar dust and its temperature. The next 
section describes the observations and the results alongwith the 
discussions are presented in section 3. 



2. Observations 

2.1. Balloon-borne observations 

The Galactic star forming region associated with IRAS 09002- 
4732 was mapped using a new 12 channel two band far infrared 
(FIR) photometer system at the Cassegrain focus of the TIFR 
100 cm (f/8) balloon borne telescope. Details of the 100 cm 




Fig. 1. The relative responses of the two far infrared bands of 
the photometer used onboard the 1 -meter balloon borne TIFR 
telescope. 



telescope system and the observational procedure have been 
described elsewhere (Ghosh et al. 1988;Mookerjeaetal. 1999). 
This 12 channel photometer was flown for the first time on 18 
November, 1993 from the TIFR Balloon Facility, Hyderabad 
(Latitude = 17?47 N , Longitude = 78?57 E). A pair of six el- 
ement (2x3 close packed configuration) composite Silicon 
bolometer arrays was used as FIR detectors and were cooled 
to 0.3 K using a closed cycle 3 He refrigerator (Verma et al. 



1993). The same region of the sky was viewed simultaneously 



in two FIR bands with near identical fields of view of 1 '.6 per 
bolometer, thus instantaneously covering an area of 6.'0 x 3. '4 
in each band. A series of filters and a Csl beam splitter define 
the two FIR wavebands. The spectral transmission of the to- 
tal system in each FIR band was determined in the laboratory 
using a Michelson interferometer and a Golay cell as the ref- 
erence detector. The resulting transmission curves for the two 
bands are shown in Fig. 1. Since the configuration of the cryo- 
genically cooled filters which define the FIR pass bands, have 
been upgraded from time to time, these bands differ from those 
in Mookerjea et al. ( |1999| ). 

The sky was chopped along the cross-elevation axis at 10 
Hz with a throw of 4.'2. The relative responses of individual 
bolometers, absolute flux calibration of the photometer as well 
as of the instrumental point spread function (PSF) were deter- 
mined from raster scans across Jupiter. The PSF also includes 
the effect of the sky chopping. The calculated effective wave- 
lengths and measured sensitivities for the two FIR bands for 
various relevant grey body spectra of different temperatures 
are presented in Table 1 . Hereafter, the two bands will be re- 
ferred to as 148 and 209 /im bands corresponding to the A e // 
for a 30 K source with A -2 emissivity. The region of the sky w 
27' x 20' , centered on IRAS 09002-4732 was mapped twice. 
The mapped region also included the IRAS Point Source Cat- 
alog (hereafter PSC) sources IRAS 08583-4719, 08589-4714 
and 09014-4736. 
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Table 1. Characteristics of the two band FIR photometer. 
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We carried out phase sensitive detection of individual 
bolometer signal using an on-board signal processing system. 
The samples from different bolometers were combined dur- 
ing the off-line data processing pipeline which reconstructs 
the telescope aspect. This uses the relative responses of the 
bolometers and a geometric model of their locations in the tele- 
scope focal plane. The resulting final chopped FIR signals were 
gridded into a two dimensional sky matrix (elevation x cross 
elevation) with a cell size of 0.'3x 0.'3. The observed signal 
map was deconvolved using an indigenously developed proce- 
dure based on the Maxim um E ntropy Method ( MEM) similar 
to that of Gull & Daniell ( |1978| ) (see Ghosh et al. |l988j for de- 
tails). The determination of the absolute aspect of the telescope 
was improved by using a focal plane optical photometer which 
detects stars (in an offset field) while the telescope scans the 
programme source. An absolute positional accuracy of ~ Q'.5 
has been achieved in the FIR maps using this method. From 
the deconvolved maps of Jupiter, the FWHM sizes are found to 
be If0xlf4and l'0xl'3 in the 148 and 209 /im bands respec- 
tively. 

2.2. Data from IRAS and ISO missions 

The data from the IRAS survey in the four bands (12, 25, 60 
and 100 /im) for the region of the sk y mapped by us were 
HIRES processed (Aumann et al. |l990| ) at the Infrared Process- 
ing and Analysis Center (IPAC[j Caltech). In the present study, 
these maps have been used for extracting sources and quanti- 
fying interband positional associations and flux densities. They 
have also been used for generating maps of dust colour tem- 
perature and optical depth. Since both IRAS 09002-4732 and 
09014-4736 appear in the IRAS Low Resolution Spectrometer 
catalog, their spectra in the wavelength range 8-22 /im have 
also been used to construct the respective SEDs. Similarly be- 
tween 45-195 /im, the data from the Long Wavelength Spec- 
trometer onboard ISO have also been used. 

3. Results and Discussion 



1 IPAC is funded by NASA as part of the part of the IRAS extended 
mission program under contract to JPL. 



3.1. 148 /im and 209 fim bands 



The MEM deconvolved TIFR maps at 148 and 209 fim of the 
Galactic star forming region associated with IRAS 09002-4732 
are presented in Fig. 2 (a & b). The insets show the decon- 
volved intensity contours for the point-like source (Jupiter), to 
quantify the achieved angular resolution. The dotted lines in 
these figures mark the boundary of the region covered by the 
telescope bore-sight. However the actual observations extend 
beyond this boundary due to the sky chopping and the effective 
size of the PSF array. As a result, the intensity maps generated 
by the MEM deconvolution extend further by 3 .'4 along in-scan 
and 6' along cross-scan directions. The limits for reliable ex- 
traction of discrete sources and structural information in the 
form of isophot contours have been established by using the 
data from two independent sets of raster scans. The dynamic 
range has been found to be ks 200 for both bands. Several re- 
gions of extended emission are seen in both the 148 and the 
209 /im maps. A total of 15 sources have been detected, de- 
tails of which are listed in Table 2 and positions marked on the 
maps. The four most prominent peaks are associated with IRAS 
PSC sources, viz., 08583-4719, 08589-4714, 09002-4732 and 
09014-4736. The listed flux densities have been obtained by 
integrating over a circle of 3' diameter. 

The positions of the global peaks in both the 148 and 209 
/im maps are identical to that of IRAS 09002-4732. The direc- 
tion of the major extended emission seen from IRAS 09002- 
4732 in both the TIFR maps [size at 10 % (50 %) of peak 
: 3.'8 x 2.'9 (1.'8 x 1.'6) at both these bands] matches quite 
well with mm-wave mole cular line maps. The CS (2-1) map 
of Zinchenko et al. ( |l995| ) shows a FWHM size of « 2' x 1 .'5 
with the major axis at almost identical position angle as the 
FIR exten sion. Higher angular resolution map by Lapinov et 
al. ( 1998 ) in the CS (7-6) line resolves the central region into 
two peaks with different LSR velocities nearly symmetrically 
displaced relative to the FIR peak. The line joining these two 
peaks also has a similar position angle. 




Fig. 2. The intensity maps for the region around IRAS 09002-4732 in TIFR bands - (a) at 148 /im with peak = 3079 Jy/sq. 
arcmin, (b) at 209 /im with peak =1611 Jy/sq. arcmin. The isophot contour levels in both (a) and (b) are 80, 40, 20, 10, 5, 2, 1, 
& .5 % of the respective peaks. The crosses mark the positions of sources detected in respective bands (see Sect. 3.1). The insets 
show deconvolved images of Jupiter in the respective bands. The contours are 60, 20, 5 & 1% of respective peaks, aligned to the 
instrumental axes for meaningful comparison. 



3.2. IRAS bands 

The HIRES processed maps in all the four IRAS bands for 
the region of present study are shown in Fig. 3 (a,b,c & d). 
The dynamic range of these HIRES maps is larger than that 
of the TIFR maps, and the isophots have been shown down to 

0. 1% of respective peak intensities. Due to this larger dynamic 
range, the diffuse emission is more prominently visible in these 
HIRES maps. 

The angular resolutions achieved in these maps (in-scan x 
cross-scan) are 0f52 x0'.63, 0f58x0.'70, l.'15xl.'35 andl.'87x 

1. '97 at 12, 25, 60 and 100 /im respectively. 

A new algorithm has been used to extract sources in the 
HIRES processed IRAS survey maps. The scheme works on the 
local maxima in the intensity maps as candidates and models 
the local background (encompassing the source of an assumed 
size, which is a variable : 1—4' ) and subtracts the same. The 
flux densities (3' diameter) for these HIRES sources are also 
presented in Table 2. 

3.3. Discrete sources 

Out of the total 15 FIR sources that have been detected from 
the TIFR maps, 12 are common to both the 209 /im and 148 
^m bands (Table 2). All these 15 sources are found to be asso- 
ciated with sources extracted from the HIRES maps in at least 
one band (13 in at least 2 HIRES bands). The positional match 



demanded is within 1 .'5 for the purpose of association. Four of 
these are associated with IRAS PSC sources whose SEDs are 
studied in greater detail (Sect. 3.6). The IRAS PSC flux densi- 
ties for these four sources are also listed in Table 2. The dust 
temperatures in the FIR, Tfir, have been computed from the 
flux densities in the TIFR bands, assuming an emissivity law 
of e\ oc A~ 2 . These are also listed in Table 2. 

An attempt has been made to determine 7, the slope of 
the IMF (dN/d(logM) oc Af 7 ), for the region around IRAS 
09002-4732, under reasonable assumptions. We assume all the 
15 sources to be nearly at the same distance based on the find- 
ing of WB89 (3 out of the 4 IRAS PSC sources in this region 
for which CO emission has been found, are at distances within 
± 10 %). If each source represents one embedded ZAMS star, 
and th eir lum inosities are translated to stellar masses using Tout 
et al. ( 1996 ), then the local slope of the IMF, 7 has been found 
to be — 1.25±g for the mass range 4-16 M@. The effect of 
incompleteness near the sensitivity limit has been avoided by 
appropriate choice of the mass bins based on the source num- 
ber distribution. This slope is similar to those found for Galac- 
tic OB associations (Massey et al. 1995) and clusters (Phelps 
& Janes |1993|). 



3.4. Diffuse emission 

Using the T( 148/209) temperature distribution for the entire re- 
gion mapped (Sect. 3.5), and the intensity maps in the same 
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Fig. 3. The HIRES processed IRAS maps for a similar region as shown in Fig. 2, in the four bands - (a) at 12 /im with peak = 
374 Jy/sq. arcmin, (b) at 25 /im with peak = 5470 Jy/sq. arcmin. (c) at 60 /im with peak = 12300 Jy/sq. arcmin. (d) at 100 /im 
with peak = 5830 Jy/sq. arcmin. The isophot contour levels in (a) & (b) are 90, 30, 10, 3, 1, .3, & .1 % and in (c) & (d) are 80, 
40, 20, 10, 5, 2.5, 1, .5, .25 & .1 % of the respective peaks. 



region, the total far-infrared luminosity has been estimated, as- 
suming the dust emissivity law to be e\ oc A~ 2 . For the re- 
gions of the intensity maps with S/N < 3 in either band, a fixed 
value of 10 K has been assumed for the dust temperature. The 
contribution of diffuse emission in this complex is found to be 
~2.8x 10 4 Lq, which is 35 % of the total emission. This does 
not include contribution from any diffuse (structureless at < 
10') component which may be associated with the Vela Molec- 



ular Ridge clouds along the line of sight. This fraction is in 
agreement with the values found for young Galactic star form- 
ing regions, viz., 25 % for Orion A region (Mookerjea et al. 
2000% 35 % for W 31 complex (Ghosh et al. [1989}. 
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Table 2. Position and flux density details of the detected sources. 
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a From HIRES processed maps unless specified otherwise. The flux densities are integrated over a circle of 3'.0 diameter. 
6 Determined using the flux densities in TIFR bands and assuming a grey body spectrum with emissivity oc A~ 2 . 
c From IRAS Point Source Catalog. 



3.5. Distribution of dust temperature and optical depth 

Although the HIRES processed IRAS maps have a much higher 
dynamic range, the angular resolution of the TIFR maps is su- 
perior to the IRAS maps at 60 and 100 /im, due to the smaller 
and circular beams employed. Taking advantage of the nearly 
identical beams of the TIFR bands at 148 and 209 /im and the 
simultaneity of observations, we have generated reliable maps 
of dust temperature and optical depth (at 200 /im, T200; trans- 
lated from 209 /im), assuming a dust emissivity of e\ oc A -2 
foll owing the sa me procedure as described in Mookerjea et 
al. ( 1999 , 2000b ). The dust temperatures are calculated based 
on the assumption that the dust is optically thin which is rea- 
sonable at these wavelengths. Considering the uncertainties in 
the flux densities measured by us and the wavelength regime 
probed by us, the temperature computed has an uncertainty of 
~ ±2K, between 15 K and 35 K which gradually increases to 
±5K at 50 K. 

Two regions with relatively stronger FIR emission have 
been selected out of the entire mapped area for presenting the 
T( 148/209) and T200 distribution. These two regions cover the 
target source IRAS 09002-4732 (hereafter, A) and the neigh- 
bouring source IRAS 09014-4736 (B). The spatial distribution 
of T(148/209) as isotherm contours and contours of constant 
T200! for the regions A and B are presented in Figs. 4 and 5 
respectively. 

The HIRES maps at 60 and 100 pm trace the distribution 
of somewhat warmer dust (30-50 K). We have used these to 
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Fig. 4. The distribution of dust temperature T( 148/209), and 
optical depth at 200 /im, T200, for the region around IRAS 
09002-4732 (A) assuming a dust emissivity law of oc A~ 2 . 
The isotherms correspond to 15 K to 50 K in steps of 5 K. 
Temperature values for a few contours are displayed to avoid 
ambiguity. The T200 contours represent 90, 70, 50, 30, 20, 10, 
5, 2 & 1 % of the peak value of 0.051. The cross marks the 
position of IRAS 09002-4732. 



generate maps of the colour temperature, T(60/100), and op- 
tical depth (noo) for the region A (around IRAS 09002-4732). 
The intensity maps at 60 and 100 /im were spatially averaged 
before computing T(60/100) and rioo in a manner similar to 



that described by Ghosh et al. (1993) for an emissivity law 
of e\ oc A" 1 . This emissivity law is more appropriate for the 
60-100 /im wavelength range. The isotherms of T(60/100) and 
contours of constant T100 are presented in Fig. 6 (a & b). 
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Fig. 5. Same as in Fig. 4, but for the region around IRAS 
09014-4736 (B). The isotherms refer to the same temperatures 
as in Fig. 4. The T200 contours represent 10, 5, 2 & 1 % of the 
global peak value of 0.05 1 . The cross marks the position of the 
TIFR source (S15) associated with IRAS 09014-4736. 



The position of the intensity peaks in both the TIFR as 
well as all four IRAS bands for the region A are identical to 
that of IRAS PSC coordinates of IRAS 09002-4732, which 
has been marked on Figs. 4 and 6. It may be noted here that 
the positions of the associated radio continuum source at 843 
MHz (Whiteoak 1992) and H2O maser source are within 15" of 



IRAS 09002-4732 (Lapinov et al. |1998|). 

The T( 148/209) temperature distribution is very asymmet- 
rical around IRAS 09002-4732. Hotter dust and a steeper tem- 
perature gradient are seen towards the north-west (NW) of 
IRAS 09002-4732 and a relatively smoother temperature dis- 
tribution in the range of 20 - 30 K is observed towards the 
south-east. The band of hotter dust is situated ~ 2'.2 to the NW 
of the peak position of the FIR intensity. The dust temperature 
drops from a value of ~ 50 K to 20 K within a distance of ss 0.8 
pc. Similar large variations in dust temperature has been seen 
for the regions associated wit h IRA S 00494+5617 and IRAS 
05327-0457 (Mookerjea et al. |2000a| ). 

The T(60/100) map shows one major peak, coinciding with 
the position of IRAS 09002-4732 and the temperature drops (to 
25 K) on all sides almost symmetrically except for the region 
to the north. The values of T(60/100) would come closer to 
T( 148/209), if the dust emissivity exponent is chosen to be - 
2 (see Fig. 6 caption). The lack of complex structure in this 
T(60/100) map vis-a-vis T(148/209) map may be due to the 
poorer angular resolution of the IRAS maps. 

The distribution of optical depth (T200) around IRAS 
09002-4732, somewhat resembles the intensity distribution : 
the peak position matches with that of Ji4g or i209- The peak 
T200 value is 0.051. There is a secondary peak in this map 
shifted towards the NW (beyond the band of higher tem- 
perature), positionally matching with relatively cooler dust 
[T(148/209)~20-25 K]. 

The T100 distribution is similar to that of T200 only near the 
peak of the intensity maps. This may not be surprising, since 
the IRAS and TIFR maps are probing dust column densities 
at different temperatures. The highest value of T100 occurs ~ 
2' south-west of IRAS 09002-4732. The value of t wo near the 
peak of the intensity maps, is 0.023. 
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Fig. 6. The (a) dust temperature T(60/100) and (b) optical 
depth at 100 /im (tioo) distribution from the HIRES maps for 
the region around IRAS 09002-4732 (A) assuming a dust emis- 
sivity law of e\ oc A -1 . The isotherms displayed in (a) corre- 
spond to the temperatures 70, 60, 50, 40, 30 & 25 K. The in- 
nermost contour is for the hottest dust (70 K). (The isotherms 
would represent 50, 45, 39, 33, 26 & 22 K if the emissivity 
index was -2). The contours in the optical depth map (b) rep- 
resent 50, 40, 30, 20 & 10 % of the global peak value of 0.1 1. 
The cross marks the position of IRAS 09002-4732. 



The T(148/209) temperature distribution in IRAS 09014- 
4736 also shows asymmetric structures (Fig. 5a). A ridge of hot 
dust runs approximately east-west ~ 1 '.5 away (south) from the 
peak of FIR intensity distribution. In addition, a local hotspot 
is seen due north-west (1.'2) from the intensity peak. The T200 
distribution shows one major peak positionally corresponding 
to ~ 20 K dust, and an isolated minor peak at the position of 
the hotspot. 

3.6. Radiative transfer modelling 

With the aim of extracting important physical parameters, an 
attempt has been made to construct radiative transfer mod- 
els for the four IRAS PSC sources IRAS 09002-4732, 09014- 
4736, 08583-4719 and 08589-4714, which have been covered 
in our FIR maps. The self consistent scheme developed by 
Mookerjea & Ghosh (1999) has been used for this purpose 
which is briefly described below. This scheme models spher- 
ically symmetric cloud with centrally embedded young stellar 
object/(s), immersed in an average interstellar radiation field. 
Although the gas-dust coupling is not considered explicitly, the 
role of the gas in the energetics of the cloud has been included 
in a self consistent manner. Two commonly used types of inter- 



stellar dust have been explored (Draine & Lee 1984, hereafter 
DL; Mezger et al. |1982[ hereafter M MP), w ith a power law 



size distribution given by Mathis et al. ( |1977[ ). For each type of 
dust, the relative proportion of graphite and astronomical sili- 
cate grains is selectable. The observational constraints include 
the SED due to the dust, angular sizes at different wavelengths 
and the radio continuum emission. 

The following parameters are explored in order to get an 
acceptable fit to all the data : (i) the nature of the embedded 
source, which could either be a single ZAMS star or a cluster 
of ZAMS stars consistent with the Salpeter Initial Mass Func- 
tion; (ii) radial density distribution law (only three power laws 
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have been explored, viz., n(r) oc r a , with a = 0, —1 or —2); 

(iii) the relative abundances of the two constituent grain types; 

(iv) total radial optical depth due to the dust (inclusive of all 
constituents) at a selected wavelength (tioq at 100 ^m); (v) the 
gas to dust ratio by mass (the predicted radio continuum emis- 
sion is sensitive to this); (vi) geometric details of the dust cloud 
(e.g. cavity size, outer size of the cloud). 

The distances to all the four IRAS PSC sources in the 
present study, have been taken to be 1.4 kpc. 

3.6.1. IRAS 09002-4732 

The observed SED of IRAS 09002-4732 has been constructed 
using TIFR, IRAS-HIRES, IRAS-LRS and ISO-LWS data 
(Fig. 7a). In ad dition , 1.3 mm continuum measurement of 
Reipurth et al. ( 1996) and near infrared data from Lenzen 



Table 4. Comparison of predicted sizes from the best fit model 
for IRAS 09002-4732 with observed FWHMs. 



( 1991 ) have been used. The total radio continuum emission has 



been found to be 3.8 Jy at 5 GHz (Caswell & H aynes |1987| ) 
from this source in a wide beam (4'4). Whiteoak (1992) found 



the flux density at 843 MHz to be 2.1 Jy. Recent high angu- 
lar resolution (synthesized beam w 2" ) measurement for the 



central region by Walsh et al. (1998) at 8.64 GHz, gives apeak 
intensity of 0.5 Jy/beam. The total luminosity and the rate of 
Lyman continuum photon generation has been estimated to be 
6.9 xlO 4 L and 7.1 x 10 47 s _1 respectively by Simpson & 
Rubin ( 199C| ). By integrating the observed SED, we estimate 
the total luminosity to be 1 x 10 5 L Q . In a study of a sample of 
ultracompact H II regions which includes IRAS 09002-4732, 



Walsh et al. ( |1999[ ) have fitted the infrared part of the SED by 
representing it as a two-component black body spectrum and 
also by a radiative transfer model. Our self consistent mod- 
elling scheme includes angular sizes at different wavelengths, 
radio continuum data as well as dust composition and grain size 
distribution. 

The best radiative transfer model for IRAS 09002-4732 has 
been identified by varying the parameters and comparing the 
following predictions with observations : (i) the SED originat- 
ing from the dust, (ii) the angular sizes at different wavebands, 
and (iii) the radio continuum emission from the gas compo- 
nent. The predicted angular sizes have been convolved with the 
relevant telescope beams, prior to comparison with observed 
FWHM sizes. The embedded source is found to be a single 
ZAMS 07 star with N Lyc = 3.6 x 10 48 s" 1 and T eff = 38,500 
K. The preferred type of dust is DL and the density distribu- 
tion is uniform. The gas to dust ratio is 100 : 1 by mass. The 
other parameters of the best fit model for IRAS 09002-4732 are 
presented in Table 3. 

The predicted SED from this model is compared with the 
observations in Fig. 7 (a). The model fits the data quite well. 
The predicted total radio continuum emissions (1.5 Jy at 843 
MHz; 3.7 Jy at 5 GHz & 3.5 Jy at 8.6 GHz) are also in rea- 
sonable agreement with observations. For our model, the H II 
region extends into the dust cloud beyond the dust-free cen- 
tral cavity (thii > Rmin)t ar, d lts s i ze 1S of sim ilar o rder 
as inferred from the 8.6 GHz map by Walsh et al. ( |l998| ). A 
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wavelengths with the measurements is presented in Table 4. 
The predicted angular sizes increase as a function of the wave- 
length and the values are in general agreement within errors. 

The IRAS LRS catalog has classified IRAS 09002-4732 
to be of type 81 based on the strengths of emission of the 
unidentified 11.3 /im feature and the [Ne II] line at 12.8 /im. 
Our inference of an embedded 07 star from the above mod- 
elling is consistent with the detection of the [Ne II] line. How- 
ever, the detection of 11.3 /im feature, generally attributed to 
be arising from single photon absorption processes of PAH 
molecules, indicates importance of non-equilibrium heating 
processes, which is beyond the scope of the present study. 

The total mass of our model cloud is 1400 M® (gas to 
dust ratio = 100) which can be compared with inferences from 
molecular line observations. The virial mass estimat e from the 
CO data of WB89 and using the relation of Evans ( 1999 ; see 
equation 8), we get the cloud mass to be 1710 Mq which is rea- 
sonably close to the mass from our model considering the un- 
certainties. On the other hand, the mass estimate from CS line 



comparison of the predicted angular sizes at mid to far infrared 



measurements is 607 Mq (Zinchenko et al. 1995). The (L/M) 
ratio of the total IR luminosity to cloud mass from the model, 
the virial mass based on CO data, and the CS core mass, are 
71, 58 and 165 L Q /M Q respectively, for IRAS 09002-4732. 

3.6.2. IRAS 09014-4736 

The FIR emission from IRAS 09014-4736 at 148 and 209 /xm is 
quite extended as evident from Fig. 2 and it is listed in the IRAS 
Small Scale Structure Catalog (SSSC; X0901-476) due to its 
extended emission in 25 and 60 /im bands. It is gratifying to 
note that the flux densities listed in the SSSC (109 Jy at 25 /im 
and 794 Jy at 60 /im) match extremely well with flux densities 
(100 and 795 Jy respectively; see S15 in Table 2) determined 
using our source extraction scheme from the HIRES processed 
maps (Sect. 3.2). The angular sizes estimated in SSSC are : 1' at 
25 /im and 1 .'6 at 60 /mi. The FWHM sizes (major x minor) 
in the 148 and 209 /im bands are 3 .'5 x 2.'7 and 3f8 x 2.'5 
respectively. The size increases with wavelength as expected. 

The observed SED for IRAS 09014-4736 has been com- 
piled using the TIFR, IRAS-HIRES and IRAS-LRS data (Fig. 
7b). No radio continuum data for this source are available in 
the literature. However, this region has been covere d by th e 843 
MHz survey with 43" angular resolution (Whiteoak 1992 ). The 
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Table 3. Parameters for the best fit radiative transfer models. 



Source a R ma x Rmin thii Tioo Luminosity Dust composition" M dust 

(pc) (pc) (pc) (1O 3 L ) Silicate:Graphite (M©) 

IRAS 08583-4719 O0 03 5.0 x 10" b 2.9 x 10~ a O02 AS) 81 : 19 0.50 

IRAS 08589-4714 0.0 0.2 8.0 x 10" 5 1.2 x 10" 3 0.05 2.4 81 : 19 0.55 

IRAS 09002-4732 0.0 1.0 3.0 x 10~ 2 5.5 x 10~ 2 0.05 100 72:28 14 

IRAS 09014-4736 -1.0 0.8 8.0 x 10~ 5 2.9 x 10~ 5 0.03 6.5 81:19 0.80 



"Dust is of DL type for all the four sources. 



quoted completeness limit of 50 mJy by Whiteoak (1992) has 
been used here as an upper limit. Modelling the SED of IRAS 
09014-4736 leads to a better fit for an density distribution 
law than for an r° law. The embedded source is found to be a 
ZAMS star of Bl type (T eff = 22,600 K). The preferred dust 
is of DL type with the radial optical depth, tioo = 0.03. The 
other parameters corresponding to the best fit model for IRAS 
09014-4736 are listed in Table 3. The spectral fit of this model 
to the observed SED is displayed in Fig. 7 (b) which is rea- 
sonably good. There is hardly any detectable radio continuum 
emission as expected for a B 1 type star. 

Near and mid-IR spectrophotometry of IRAS 09014-4736 
shows re combi nation line emiss ion of hydr ogen ( Persson & 
Campbell [l987| ; Beck et al. |l99l| ; Porter et al. |1998| ). The LRS 
spectrum of this source has been classified as type 80 based 
on the detection of 11.3 /im feature and absence of any fine 
structure line from ionized heavier elements, typical of low ex- 
citation H II regio n. Det ection of several addit ional PAH fea- 
tures (Cohen et al. 1989; de Muizon et al. 199C; Zavagno et al. 



1992) further supports that IRAS 09014-4736 has an embedded 



YSO. Although IRAS 09014-4736 is classified as an unidenti- 
fied source in the IRAS PSC, it is situated within sa 20" of the 
small optical nebula, BB W 225 (size < 2' ), identified by Brand 
et al. ( |l986| ). The embedded YSO (Bl) inferred from our mod- 
elling is consistent with all of the above. However, the actual 
geometry of the source has to be far from spherically symmet- 
ric, in order to reconcile with the size of the optical nebula and 
the near infrared data. Similar conclusions have been arrived at 



by Beck et al. ( |1991| ) for IRAS 09014-4736. 

The only available molecular line observation from this 
source is the CO emission (WB89). We estimate the molec- 
ular gas density and the virial mass for the cloud associated 
with IRAS 09014-4736 to be 5.2 x 10 3 cm" 3 and 395 Mq re- 
spectively. The mass of the model cloud is 80 Mq for a gas to 
dust mass ratio of 100. The latter is not constrained well due to 
the lack of radio continuum emission. The (L/M) ratio corre- 
sponding to the model and the virial mass of this cloud are 8 1 
and 16 Lq/Mq respectively. 

3.6.3. IRAS 08589-4714 

In addition to the meas ureme nts presented here, the near-IR 
data of Lorenzetti et al. ( 1993 ) have been used to construct the 
SED of IRAS 08589-4714 (Fig. 7c). The parameters of the best 



fit radiative transfer model are presented in Table 3. A constant 
density cloud with DL type dust grains is implied. The total 
radial tioo = 0.05 and the model cloud mass is 55 Mq. Since 
the energizing stellar source is detectable in the near-IR, this is 
an example of very evolved stage of star formation, where only 
a small fraction of the original interstellar matter is left around 
the star at present. 

The source IRAS 08589-47 14 has been detected in CO and 



CS sur veys by Wouterloot & Brand ( 1989 ) and Bronfman et al. 
( 1996| ) respectively. Whereas our model gives a value for the 
L/M ratio to be « 44 Lq/Mq, the virial mass of the cloud 
estimated from the CO data of WB89 indicates a value of 3.7 

Lq/Mq. 

3.6.4. IRAS 08583-4719 

IRAS 08583-4719 is contained in the catalog of bright rimmed 
clouds which are candidates for star formation by radiation 
driven implosion (Sugitani & Ogura 1994). Unfortunately no 
additional observational data pertaining to a SED are avail- 
able other than those presented here. The physical size estimate 
has been taken from Sugitani & Ogura (1994). Once again a 
constant density dust (DL) cloud gives the best fit to the SED 
(Fig. 7d). However, the fit to the observations even for this best 
model is not very satisfactory. The parameters for this model 
are listed in Table 3. The mass of the model cloud is 50 Mq 
leading to a (L /M) ratio of 80 L Q /Mq. 



3.6.5. Inferences from the modelling 

The DL type of dust is favoured for all the four sources mod- 
elled here. It may be noted that the emissivity index in the 
100 - 200 fim range for DL type dust (for both silicate and 
graphite), is w -2, consistent with our choice for generation of 
T(148/209) maps (Sect. 3.5). Three of these four sources sug- 
gest a constant density dust envelopes around the YSOs. This 
conclusion seems to be very common for several Galactic star 
forming regions : e.g. for IRAS 18 314-0 720, 18355-0532 and 
18316-0602 (Mookerjea & Ghosh, |l99§); IRAS 00338 +6312, 
04004-5 1 14 and RAFGL 5111 (Mookerjea et al. |l999|); IRAS 



00494+5617 and 05327-0457 (Mookerjea et al. 2000a). The 



following comments can be made by comparing the parame- 
ters for the four sources modelled here : The dust cloud corre- 
sponding to IRAS 09002-4732 has the largest size of the central 
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Fig. 7. Comparison of the predicted spectral energy distributions (SEDs) from the best fit models of (a) IRAS 09002-4732, (b) 
IRAS 09014-4736 (c) IRAS 08583-4719 & (d) IRAS 08589-4714 with observations. The filled circles represent the TIFR and 
IRAS-HIRES data. The dashed lines are IRAS-LRS spectra (for IRAS 09002-4732 and IRAS 09014-47 36) an d the dash dot line 
is the ISO-LWS spectrum of IRAS 9002- 4732. The triangles denote 1.3 mm data from Reip urth e t al. ( |1996[ ) for IRAS 09002- 
4732 and near-I R data from Lenzen ( |1991| ) for IRAS 09002-4732, from Persson & Campbell ( |1987| ) for IRAS 09014-4736, from 
Lorenzetti et al. ( |1993[ ) for IRAS 08589-4714. The solid line denotes the prediction of the best fit model. See text and Table 3 for 
details of the model parameters. 



cavity, which contains the most luminous embedded source. 
Except IRAS 09014-4736, for all other three sources, the Hll 
regions extend well into the respective dust clouds. This em- 
phasizes the role of dust grains on the emission of recombi- 



nation lines and radio continuum. The radial optical depth due 
to the dust ranges between 0.02 to 0.05 which is slightly on 
the lower side compared to IRAS 00494+5 617 (M ookerjea et 
al. EoOOal), 00338+6312 (Mookerjea et al. (19991). The inter- 
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stellar dust composition is found to be nearly identical in all 
the four sources modelled, viz., Silicate dominated. The ratio 
between the graphite and the silicate grains, in general vary 
drastically among the Galactic star forming regions (Mooker- 



jea et al. 1999, 2000a). Hence, nearly identical dust composi- 



tion found here further supports the hypothesis that the entire 
complex around IRAS 09002-4732 mapped by us is physically 
associated. 

The L /Al ratio for the four sources studied here, as found 
from our radiative transfer modelling, falls in a narrow range 
of 44 - 81 Lq/Mq. For three of these for which CO mea- 
surements are available, this ratio based on virial masses, fall 
in somewhat broader range (3.7 - 58 Lq/Mq). These can be 
compared with the results from other studies : (i) A study sim- 
ilar to the present one, for the W 31 complex has lead to a 
val ue of (L/M) in the range of 1.6 - 3.6 Lq/Mq (Ghosh et 
al. |1989| ); (ii) the average value for a large sample of clouds 
using virial mass based on CO measurements, is ~ 4 Lq/Mq 
with a spread exceeding a factor of 100 (Mooney & Solomon 



1988 , Evans 1991); (iii) for denser cores studied in CS emis- 
sion (Plume et al. 1997 ) give ^190 Lq/Mq with a spread of 
a factor of 15. Considering the uncertainties arising due to the 
fact that the different tracers are sampling different components 
of the ISM (at possibly different spatial scales), the star form- 
ing complex associated with IRAS 09002-4732 can be termed 
as typical of the Galaxy. 

4. Summary 

The Galactic star forming region associated with the IRAS 
source 09002-4732 and its neighbouring region (w 0.15 sq. 
deg.) has been mapped in two far infrared bands at 148 and 209 
fim with an angular resolution of ~ 1'. A total of 15 sources 
have been detected including 4 IRAS PSC sources. Several 
of these are well resolved in both these bands. Diffuse emis- 
sion contributes 35 % of the total FIR emission at 209 /im. 
Taking advantage of simultaneous observations using identi- 
cal beams in these two bands, reliable maps of dust colour 
temperature [T( 148/209)] and optical depth (T200) have been 
generated which show many structures. Under reasonable as- 
sumptions, the IMF slope for this region has been found to be 

-1.25±°:65 for the mass ran 8 e of 4-16 Mq. 

The HIRES processed IRAS survey data for the same re- 
gion at 12, 25, 60 and 100 /im have also been used for compari- 
son as well as extracting discrete sources positionally matching 
with those detected at 148 and 209 /im. The spatial distribution 
of T(60/100) and tioo have also been generated for the region 
locally around IRAS 09002-4732 for comparison. 

Self consistent radiative transfer calculations has been car- 
ried out for IRAS 09002-4732, 09014-4736, 08583-4719 and 
08589-4714 through spherical dust-gas clouds. The observed 
infrared sub-mm SED, angular sizes and the radio continuum 
data have been used to constrain the parameters of the models. 
Whereas for IRAS 09014-4736 an r -1 radial density distribu- 
tion law is preferred, all the other three sources favour an r° 
law. The geometrical details of the dust-gas clouds, the dust 



compositions and optical depth, etc. have been quantified by 
these models. The luminosity per unit mass for these sources 
lie in a rather narrow range of 44-81 Lq/Mq. 
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